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ABSTRACT 


Hallen  type  integral  equations  are  formulated  for  an  arbitrarily 
bent  monopole  above  an  infinite  perfect  ground  plane  and,  by  numerical 
solution  techniques,  current  distributions  are  obtained.   The  coupled 
integral  equations  are  founded  upon  the  familiar  electric  scalar  poten- 
tial and  magnetic  vector  potential.   Numerical  results  are  presented 
for  several  cases  including  different  bent  angles  and  are  compared  where 
possible  with  available  data. 
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I.   INTRODUCTION 

A  linear  antenna  is  one  in  which  a  "long  dimension"  can  be  discerned; 
for  example,  a  long  circular  cylinder  the  length  of  which  is  much  greater 
than  the  diameter  is  looked  upon  as  a  linear  antenna.  Initially,  the  em- 
phasis will  be  limited  to  the  so-called  cylindrical  antenna  driven  at  its 
mid-point  by  a  voltage  source  V  applied  to  a  vanishingly  short  gap  at  the 
antenna  center.  The  antenna  will  be  considered  to  reside  in  a  lossless 
homogeneous  and  isotropic  medium  of  infinite  extent. 

A  thin  wire  antenna  is  commonly  defined  as  a  linear  antenna  which  has 
a  radius  much  less  than  a  wavelength.   There  are  several  approximations 
which  may  be  employed  in  the  analysis  of  thin  wire  antennas. 

(1)  The  current  on  the  cylindrical  surface  is  assumed  to  be 
entirely  in  the  axial  direction. 

(2)  The  current  distribution  is  assumed  to  be  constant  around 
the  periphery  of  the  wire  at  a  given  cross-section  and  hence  can  be 
approximated  as  a  filament  of  current  on  the  axis  of  the  cylinder. 

(3)  The  boundary  conditions  are  applied  only  to  those  components 
of  the  total  electric  field  which  are  tangent  to  the  wire  surface  and 
parallel  to  the  wire  axis. 

Many  investigators  have  considered  the  theoretical  problem  of  deter- 
mining the  current  induced  on  a  thin  wire  scattering  antenna.  [Ref.  1,2,3. 
These  investigations  have  been  extended  to  include  the  treatment  of  ar- 
bitrary configurations  of  thin  wire  scattering  antennas  intersecting  at 
a  single  junction.  [Ref.  4].   Three  particular  configurations  for  which 
numerical  results  have  been  obtained  are  a  set  of  nonintersecting  thin 


straight  wires,  [Ref.  5],  a  set  of  perpendicular  intersecting  thin  straight 
wires,  [Ref.  6],  and  a  pair  of  skewed  crossed  wires.  [Ref.  7]. 

One  general  approach  to  finding  the  current  on  thin  wire  antennas  is 
to  obtain  a  set  of  integral  equations  which  involve  the  unknown  current 
distributions.   Using  the  equation  of  continuity,  one  may  determine  ex- 
pressions for  the  vector  and  scalar  potentials  in  terms  of  the  current 
distribution.   It  is  in  obtaining  these  integral  equations  that  a  choice 
can  be  made. 

Following  the  Pockington  equation  formulation,  the  expressions  for 
the  vector  and  scalar  potentials  are  used  along  with  Maxwell's  equations 
to  derive  a  system  of  integro-dif ferential  equations  for  the  unknown 
current  distribution.   An  auxiliary  scalar  function  is  defined  and  intro- 
duced into  the  system  of  integro-dif ferential  equations.   The  solution  of 
these  now  ordinary  differential  equations  subject  to  appropriate  boundary 
conditions  gives  the  desired  set  of  integral  equations  for  the  unknown 
current  distribution.  , 

If  the  Hall£n  integral  equation  approach  is  used,  one  independently 
establishes  two  expressions  for  a  given  component  of  the  vector  potential 
in  the  direction  of  a  wire  element  and  these  expressions  are  evaluated  on 
the  wire  surface.   One  expression  is  obtained  from  the  solution  of  an 
inhomogeneous  differential  equation  satisfied  by  the  component  of  vector 
potential. 

Another  expression  for  vector  potential  is  obtained  directly  from  the 
familiar  vector  potential  itself.   These  two  expressions  must  be  equal 
and,  hence,  are  equated.   Next,  the  boundary  condition  that  the  tangential 
component  of  electric  field  on  the  surface  of  the  wire  element  must  be 
zero  is  applied.   Now  the  arbitrary  constants  of  integration  which  arise 


in  the  solution  of  the  inhomogeneous  differential  equation  are  evaluated. 
The  constants  of  integration  must  have  assumed  values  which  are  in  com- 
pliance with  all  other  required  conditions.  [Ref.  8]. 

There  have  been  numerous  methods  proposed  for  solving  the  integral 
equations,  however  obtained.   These  methods  have  fallen  into  a  few  broad 
catagories: 

(1)  Iteration  [Ref.  9,  10]. 

(2)  Expansion  of  the  current  in  a  Fourier  series  [Ref.  11]. 

(3)  The  variational  method  [Ref.  12,13,14]. 

(4)  The  Wiener-Hopf  techniques  for  the  case  of  the  long 
antenna  [Ref.  15]. 

(5)  Numerical  methods  [Ref  1,  16]. 

One  effective  numerical  method  of  solving  the  integral  equation  is 
the  so-called  direct  integration  technique  (Method  of  Moments).  [Ref.  16, 
17].   Basically  the  technique  utilizes  an  expansion  of  the  unknown  func- 
tion in  some  convenient  manner.   This  expression  may  be  a  finite  sum  of 
piecewise  continuous  functions  which  are  non-zero  over  some  small  inter- 
val of  definition  and  vanish  elsewhere.   A  substitution  of  this  expres- 
sion for  the  unknown  function  into  the  system  of  integral  equations 
reduces  them  to  a  system  of  linear  algebraic  equations  for  the  unknown 
coefficients.   These  linear  algebraic  equations  may  be  solved  by  conven- 
tional techniques. 

The  purpose  of  this  thesis  is  to  formulate  the  coupled  integral 
equations  for  an  arbitrarily  bent  monopole,  above  an  infinite,  perfectly 
conducting  ground  plane,  in  a  Hallen-type  integral  which  is  founded 
entirely  upon  the  familiar  scalar  potential  and  magnetic  vector  potential, 
These  integral  equations  are  then  solved  using  the  direct  integration 


technique.   Throughout  the  thesis  it  is  found  that  the  coupled  integral 
equations  obtained  by  the  Hallen  approach  can  be  readily  obtained  with- 
out employing  an  auxiliary  scalar  potential.   How  the  various  arbitrary 
constants,  which  arise  in  the  development,  are  related  to  the  electric 
scalar  potential  is  clearly  shown. 
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II.   DERIVATION  OF  MATRIX  EQUATION 

In  this  section  the  system  of  integral  equations  for  the  current 
distribution  on  an  arbitrarily  bent  monopole,  above  an  infinite,  perfect- 
ly conducting  ground  plane,  driven  at  its  base  by  a  delta  gap  generator 
is  derived.   From  the  equations,  an  operator  is  defined  such  that  the 
method  of  moments  can  be  applied  to  find  an  approximate  representation 
of  the  operator  by  a  matrix.   The  solution  for  the  desired  current  dis- 
tribution is  obtained  by  inverting  the  matrix. 

A.   FORMULATION  OF  THE  PROBLEM 

Figure  1  defines  the  coordinate  system.   The  antenna  consists  of  two 
thin  wire  segments  with  a  single  junction.   The  segments  have  radii  a 
and  a  and  lengths  h   and  h  .   The  primary  source  is  the  field  in  the 
gap  which  is  maintained  at  V  /A  volts  per  meter  by  the  generator.   The 
induced  source  is  the  current  on  the  dipole.   The  solution  is  carried 
out  assuming  the  wires  are  perfect  conductors.   With  this  assumption  the 
tangential  component  of  the  total  electric  field  at  the  surface  of  each 
wire  segment  must  vanish.   Due  to  axial  symmetry  of  each  wire  segment,  the 
phi-component  of  each  electric  field  must  vanish.   Hence  for  each  wire 
segment  the  following  is  obtained: 

E^r  )  +  EP(7  )  =  0  (1) 

z   s     z   s 

EX> +  EsQ  ■  °  (2> 

where  E  (r  )  is  the  tangential  component  of  the  induced  electric  field 

at  the  surface  of  the  z-directed  wire  and  EP(r  )  is  the  tangential 

z   s 
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Figure  1.   Antenna  Geometry 
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component  of  the  primary  electric  field  at  the  surface  of  the  z-directed 
wire.   Equation  (2)  is  similarly  obtained  for  the  s-directed  wire. 

B.   DERIVATION  OF  THE  COUPLED  INTEGRAL  EQUATIONS 

The  Halle'n  equations  are  derived  from  a  classical  viewpoint  for  the 
case  of  a  cylindrical  antenna  driven  at  its  center  by  a  so-called  slice 
generator.  [Ref.  9,  18,  19].   The  general  procedure  for  developing  a 
Hallen  type  equation  is  as  follows: 

(1)  Determine  the  magnetic  vector  potential  component  from  the 
potential  integral  equation  for  A. 

(2)  Determine  the  magnetic  vector  potential  component  along  the 
cylinder  from  the  inhomogeneous  differential  equation  which  the  magnetic 
vector  potential  component  satisfies. 

(3)  Write  the  desired  integral  equation  by  equating  the  two  ex- 
pressions obtained  in  steps  one  and  two. 

(4)  Apply  appropriate  boundary  conditions  on  the  surface  of  the 
cylinder  to  the  equation  obtained  in  step  three. 

1.   Vector  Potentials 

Consider  the  antenna  shown  in  figure  1.   Assuming  harmonic  time 
dependence  (e   ) ,  the  vector  potential  due  to  the  antenna  current  is 


4tt 


-jk | r-r' | 
J(r')  - dv(r')  (3) 


t                 -jk|p-p'| 
A(P)  =  -f-       J(p')  e     dv(p')  (4) 

£.        |P-P'I 
Due  to  the  geometry  the  vector  potential  has  only  a  z-directed  and  s- 
directed  component  which  can  be  written  on  the  axis  of  each  wire  segment 


13 


as 


A  (z)  -  •£■ 

Z  4  71 


-JkR 


I  (z«)  ^— 
z       R 


dz» 


(5) 


A  (s)  =  -y- 
s      4tt 


-JkP 


Vs'>  ^T 


dsf 


(6) 


where  R  =  |r-r'|  and  P  =  jp-p'|.   R  and  P  were  evaluated  at  r'  and  p' 

equal  to  the  radius  of  the  respective  wire  segments  with  <f> '  equal  to  zero, 

2 .   Non-Homogeneous  Differential  Equations 

It  is  now  seen  that  A  (z)  and  A  (s)  can  not  be  determined  di- 

z         s 

rectly  from  equations  (5)  and  (6)  without  knowledge  of  I(z')  and  I(s'). 
A  differential  equation  is  now  derived  which  can  be  solved  for  A  (z)  and 
Afs) 


E  =  -juoA  + 


grad  (div  A) . 


(7) 


The  electric  fields  produced  by  the  currents  are,  according  to  equation 
(7), 

1/(7)  -  j,.Az(r)Gz  +  -^    grad  [div(^ff)U,)]  (8) 

V<F)  =  -J»As(p)Hs  *  -J+-     grad  [di^A.ff)5.)]-  <9> 

Having  applied  the  gradient  operation  in  circular  cylindrical  coordi- 
nates and  with  axial  symmetry  eliminating  the  phi-component  of  E   (r) 

and  E   (p),  the  fact  that  the  total  tangential  electric  field  on  the 

surface  of  the  antenna  must  vanish  is  applied.   The  following  results 
are  obtained: 


A    (r   )] 
_   z      s  J 


+  k2  a  (7  )  sJ!Le  p(7  )  -#■ 


z      s 


uu 


A(r) 


dz  Ldy     y     s 


(10) 


[w1 


j+   k     Vps>  = 


^  E   P(F  ) 

(JJ        s    vrs 


fA 

_ow     w 


^s>]> 


(11) 
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where  r   and  p   denote  any  point  on  the  z  and  s-directed  wires  respec- 
tively.  Equations  (10)  and  (11)  are  two  coupled  partial  differential 
equations  because  of  the  last  term  on  the  right  in  each  equation.   A 
Hallen  solution  of  these  two  coupled  equations  is  needed. 
3.   Solution  of  Coupled  Differential  Equations 
Equations  (10)  and  (11)  are  of  the  form 


dx 


(12) 


The  general  solution  to  equation  (12)  was  contained  in  Ref.  [23].   It 
was  found  to  be 


f  (x)   =   A  cos   kx  +   B    sin  kx  + 


■>x 


~x   =q 


TT/        '\      A.      dV  (X  '   )"]  .  ,      ,  >\J        ' 

U(x    )   +  — ^  N  >    f'\    sin  k(x-x    )dx 

_  OX  J 

(13) 


The  integral  portion  of  equation  (13)  is  zero  at  x  equal  q  and  its  deriv- 
ative is  zero  at  this  position.   One  integration  by  parts  results  in 


r»x 


-  v  / —  sin  k(x-x/)dx/  -  -V(q)  sin  k(x-q) 

/   /   ox 


~x  =q 


px 


+  k 


V(x')  cos  k(x-x')dx 


°x  =q 
which,  when  put  into  equation  (13),  gives 


(14) 


f  (x)  =  A  cos  kx  -f  B  sin  kx  -  —  V(q)  sin  k(x-q) 


r»x 


+ 


r»x 


U(x/)  sin  k(x-x')dx  + 


°x  =q 


°x  =q 


V(x  )  cos  k(x-x  )dx 
(15) 


Since  V(q)  is  fixed,  the  first  three  terms  are  combined  and  new  constants 
obtained  are 


a  -  A  +  —  V(q)  sin  kq,  and 
b  =  B  —  —  V(q)  cos  kq. 


(16a) 
(16b) 
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Combining  equation  (16)  with  equation  (15)  gives 


r>x 


f  (x)  =  a  cos  kx  +  b  sin  kx  +  — 


U(x  )  sin  k(x-x')clx 


°x  -q 


ox 


+ 


V(x  )  cos  (x-x  )dx  . 


(17) 


°x  =q 
With  equation  (17)  available  the  solution  to  equations  (10)  and  (11)  is 
determined  with 


-^  Ay(rs)  =  V(x);  ijp  E^Cr,)  =  U(x)    for  equation  (10) 

a  -k2        - 

and   -3-  A  (p  )  =  V(x);  1—  E  P  (p  )  =  U(x)    for  equation  (11) 

OW    w          i>  UU     b      S 


which  gives 


jk 


A  (z)  =  a   cos  kz  +  b   sin  kz  +  —      E  (n)  sin  k(z-n)dn 
z       z  z  O)  ,     z 

4i=z0 


1  Pz    T  9     n 

-r-       T~  A  (n)   cos  k(z-n)dn    and 
k  .     Loy  y   J 


(18a) 


ik 

A  (s)  -  a   cos  ks  +  b   sin  ks  +  _      E  (n)  sin  k(s-n)dn 
s       s  s  uu  I  _    s 

o 


k 


4l=S. 


5_ 
_Bw 


A  (n)   cos  k(s-n)dn. 


(18b) 


The 


c)y 


A  (z) 

-  y  ■ 


means  A  (r)  is  differentiated  with  respect  to  y  and  then 


has  to  be  evaluated  along  a  z-directed  line  on  which  it  is  dependent. 

h   r    ~i 

This  is  similarly  true  for  r—  A  (s)  .   The  arbitrary  constants  of  equa- 

ow  L  w   J 

tion  (18a)  are  a  ,  b  and  z   and  any  one  of  these  three  may  be  specified 
directly,  but  the  other  two  must  comply  with  physical  boundary  condi- 
ditions.   This  applies  to  equation  (18b)  and  constants  a  ,  b   and  s 
also. 
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4.   Scalar  Potential 

With  equations  (18a)  and  (18b)  one  faces  the  problem  of  what  to 
do  about  the  constants  of  integration.   One  can  enforce  the  boundary 
condition  that  the  scalar  potential  must  be  continuous  at  each  wire  seg- 
ment junction  to  eliminate  a  constant  or  relate  two  constants.   To  obtain 
such  an  expression  is  the  next  task. 

Differentiating  equation  (18a)  with  respect  to  z  and  (18b)  with 
respect  to  s  according  to  Leibnitz's  rule  [Ref.  24]  gives 


t-  A  (z)   =  -ka  sin  kz  +  kb   cos 
3z   z  z  z 


kz  +  —  {^     E  (n)  sin  k(z-n)dn} 

dZ  l  U  J  Z  J 


n=z 
o 
z 


^yz)  +k 


o 
and 


)]  sin  k(z-n)dn  (19a) 

y 
n=z 


L~  A  (n: 

.  dy  y 


s 


3   'k  f 
A  (s)   =  -ka  sin  ks  +  kb   cos  ks  +  -r—  {-—  E  (n)  sin  k(s-n)dnl 
s  s  s  9s  L  a)  j  s  J 


3s   s 

n=s 


o 


'ftp-  A  (n] 
J  yw  w 


-  —  Aw(s)  + 


n=s 
o 

The  Lorentz  condition  relates  the  magnetic  vector  potential  A  to  the 

electric  scalar  potential  and  yields 

2 

-  J^H>(7)   =  ~  A   (7)  +  -ri  A   (7)        and  (20a) 

to  3y     y  3z     z 

2 

-  ^>(p^   =  ~  AT  G)  +  -ri  A(p")    .  (20b) 

to  3ww  3ss 

Rearrangement  of  equations  (19a)  and  (19b)  gives 

2  z 

-  J d>  (z)  =  -ka  sin  kz  +  kb   cos  kz  +  -r—  \—       E  (n)  sin  k(z-n)dn} 

u  z  z  3zlwjz 

n=z 
z  o 

+k    [[—  A   (n)]     sin  k(z-n)dn  (21a) 

J  dy     y 

n=z 
o 
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and 


d>  (s)  =  -ka  sin  ks  +  kb  cos  ks  +  — -  {—  E  (n)  sin  k(s-n)dn) 

u  s  s         3s  l  u   s  J 


j 
n=s 


+k 


[—  A  (n)]  sin  k(z-n)dn, 

dW   W 


(21b) 


n=s 


which  is  valid  along  the  z-  and  s-directed  lines  respectively. 

(j)  (z)  is  the  total  scalar  potential  along  the  z-  directed  wire 

g 
and  <J>  (s)  is  the  total  scalar  potential  along  the  s-directed  wire,  both 

of  which  are  due  to  all  charges  on  the  wire  segments,  but  the  potentials 

are  expressed  in  terms  of  the  constants  of  integration.   The  partial 

derivatives  expressed  in  equations  (21a)  and  (21b)  must  now  be  evaluated. 

Going  back  to  equation  (13)  it  is  seen  that  the  integral  portion 

of  (13)  is  zero  at  x  equal  q  and  its  derivative  is  zero  at  x  equal  q. 

Now  in  the  derivation  these  forms  are  used 


x  x 

f  { 

t-   |  U(x')  sin  k(x-x')dx"  +   V(x')  cos  k(x-x')dx"  . 


(22) 


x   =q  x  =q 

Their   relationship   to   the   antenna  problem  is   given  below  equation    (17). 
Now  applying   these    facts    the    result    is 


1  { J]i 

3z     "■     to 


~  I2-       E    (n)    sin  k(z-n)dn|    =  0      for   all   z 


n=z 


3s    i    oj 


[^       E    (n)    sin  k(s-n)dn}    =  0      for   all  s    . 


(23a) 


(23b) 


n=s 


This    result   implies    that   one  may   choose   s      and   z      as   desired,   but   at    the 

o      o 

same  time  doing  this  in  a  manner  to  lessen  subsequent  calculations.   Once 

s   and  z   are  determined  the  other  two  constants  of  integration  in  each 
o      o 

equation  must  meet  boundary  conditions. 
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C.   BOUNDARY  CONDITIONS 

The  application  of  these  necessary  boundary  conditions  to  the  cur- 
rents on  the  thin  wires  should  provide  the  additional  equations  to  yield 
a  unique  solution  for  the  current  distributions. 

1.  At  The  Junction 

As  has  been  stated  by  Taylor  [Ref.  4],  there  are  essentially 
two  conditions  which  must  prevail  at  all  junctions: 

(1)  Kirchoff's  current  law  must  hold. 

(2)  Continuity  of  the  scalar  potential  must  exist. 

The  first  of  these  requires  that  there  be  no  buildup  or  accumulation  of 
charge  at  the  junction.   Stated  differently,  the  total  charge  entering 
the  junction  must  be  equal  to  the  total  charge  leaving  the  junction. 
The  second  condition  requires  that  the  potential  on  each  of  the  wires  at 
the  junction  be  the  same.   There  can  be  no  jump  or  discontinuity  in  sca- 
lar potential  at  a  wire  junction. 

2.  At  The  Free  Ends 

According  to  previous  considerations  of  thin  wire  theory  [Ref.  1, 
2,4],  the  current  must  vanish  at  the  free  ends  of  all  wires  in  the  con- 
figuration. 
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III.   SIMPLIFICATIONS  ARISING  FROM  GEOMETRICAL  SYMMETRY 

Any  symmetries  in  the  geometry  of  the  antenna  elements  are  reflected 

in  the  structure  of  the  generalized  impedance  matrix  Z    .   Judicious 

mn 

use  of  these  symmetries  can  result  in  considerable  reduction  in  both 
computing  time  and  the  number  of  memory  locations  required  for  a  solution, 

A.  IMAGE  SYMMETRY 

Referring  to  figure  1  the  antenna  has  image  symmetry  about  the  plane 
z  =  0  through  the  midsection  of  the  antenna.   The  boundary  condition  at 
a  perfect  electric  conductor  is  zero  tangential  components  of  electric 
field.   An  element  of  source  plus  an  "image"  element  of  source,  radiating 
in  free  space,  produces  zero  tangential  electric  field  over  the  plane 
through  the  midsection  of  the  antenna.   According  to  uniqueness  concepts, 
the  solution  to  this  problem  is  also  the  solution  for  a  current  element 
adjacent  to  a  plane  conductor.   The  necessary  orientation  and  excitation 
of  image  elements  is  summarized  in  many  texts.   The  main  point  though 
is  that  symmetry  must  be  maintained  in  the  image  problem.   Once  image 
theory  is  used  and  the  resulting  symmetry  is  recognized,  all  that  is 
needed  is  to  find  the  kernel  for  each  "image"  current  element  and  com- 
bine it  algebraically  with  its  corresponding  "non-image"  current  element 
kernel  which  reduces  the  number  of  required  equations  by  a  factor  of  two. 

B.  TOEPLITZ  PROPERTY 

The  Toeplitz  property  of  the   Z    matrix  can  be  used  to  realize 

mn 

considerable  economy  of  computing  time.   An  element  of  the  matrix  depends 
only  on  the  absolute  value  of  the  distance  between  the  source  point  m 
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and  the  field  point  n  if  there  is  no  bend  involved.   Since  there  are 
only  N  distinct  distances  involved,  only  N  different  quantities  need 
be  calculated  for  a  matrix  containing  N  x  N  elements.   These  N  quanti- 
ties are  calculated  first  and  the  computer  is  instructed  to  arrange 
them  in  the  proper  locations  to  form  the  complete  matrix. 

The  antenna  has  a  bend  so  only  certain  portions  of  the  matrix  can 
be  obtained  .using  this  property.   It  is  the  image  of  the  bend  that  pre- 
vents the  entire  matrix  from  exhibiting  this  property. 
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IV.   NUMERICAL  SOLUTION 

One  of  the  most  general  and  most  useful  techniques  for  solving  pro- 
blems in  antenna  theory  is  the  direct  integration  technique.   A  detailed 
discussion  of  the  technique  has  been  presented  by  Harrington  [Ref.  16], 
Since  this  basic  technique  is  employed  throughout  the  remainder  of  this 
thesis,  a  short  quantitative  discussion  using  the  notation  of  Harrington 
is  included  at  this  time. 

A.   METHOD  OF  MOMENTS 

The  equation  to  be  solved  is  of  the  foi"m 

L(f)  =  g  (24) 

where  L  is  a  linear  operator,  g  is  the  tangential  component  of  the  applied 
field  on  each  conductor  and  f  is  the  unknown  current  to  be  determined. 
Using  the  terminology  of  Harrington  [Ref.  17]  it  is  evident  from  the 
uniqueness  theorem  that  the  problem  is  deterministic.   Since  L  and  g  are 
given,  it  is  a  problem  of  analysis. 

A  unifying  concept  for  the  solution  of  such  problems  is  afforded 

by  the  method  of  moments.   In  order  to  apply  the  method,  it  is  necessary 

to  define: 

(1)  A  set  of  expansion  or  basis  functions  f „ ,  f „,...,  f  in  the 

1   2       n 

domain  of  L.   These  known  functions  must  be  linearly  independent  and 
chosen  such  that  a  reasonable  approximation  to  f  can  be  obtained  by  a 
finite  series  expansion  of  the  form 

N 

f  =  J  a   f  (25) 

fc,  n  n 
n=l 
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where  the  a  's  are  unknown  constants  to  be  determined.   If  equation  (25) 
n 

is  substituted  into  equation  (24),  one  obtains 


N 


I   «  L(f  )  -  §  ' 


n=l 


n   n 


(26) 


(2)  A  set  of  weighting  or  testing  functions  W..  ,  w\,...,W  ,  W  ,  ,, 

1        2  m        m+1 

. . . ,W„  in   the   range   of   L  and   linearly   independent  must   be    chosen.      If 
'    N 

the  same   set   of   functions   is    used   for  expansion   and   testing,    the   pro- 
cedure  is   known   as   Galerkin's   method  and   is   equivalent   to   the   Rayleigh- 
Ritz   variational  procedure    [Ref.    17,    25]. 

(3)  An   inner  product   suitable   for   the  problem.      For  most   electro- 
magnetic  field  problems,    a  suitable   inner  product   is 


n, 


<W,f> 


f(n)W(n)dn' 


(27) 


n. 


where  n   and  n„  are  chosen  over  the  region  of  interest.   The  resulting 
expression  is 


N 


I   a      <W  ,L(f  )>  -  <W  ,g>  ,   m  =  1,2,...,N  . 
L,  n   m    n       m'fo   '       '  '    ' 


(28) 


n=l 


Inspection  of  equation  (28)  shows  that  it  can  be  written  as  a  matrix 
equation  by  defining  the  following  matrices 


mn 


<W1,Lf1>  <W1,Lf2> 
<W2,Lf1>  <W2>Lf2> 


W  <«N.«2> 


<W2'LV 


W 


(29) 


[«nl   = 


,  and   [V  ]  = 
m 


<w1,g> 
<w2,g> 


<wN,g> 


(30) 
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The  matrix  equation  may  now  be  written  as 

l\J  K]'  [V  .  (3D 

By  analogy  with  multiport  network  theory,  [z   ]  is  called  a  generalized 

mn 

impedance  matrix  and  [V  |  a  generalized  voltage  matrix.   Only  if  [z 

m  mn 

is  nonsingular  does  its  inverse  exist  and  the  desired  solution  is 

[a  ]  -  [z  J"1  [V  ]  .  (32) 

n      mn      m 

Since  the  a  's  are  then  known,  the  solution  for  the  unknown  current  f 
n 

is  determined  by  equation  (25) . 

The  expediency  of  obtaining  a  solution  for  f  depends  not  only  on  L 

and  g,  which  are  usually  determined  by  the  physical  problem,  but  also 

on  the  functions  f   and  W  .   Extreme  care  should  be  taken  in  choosing 

n      m 

these  two  sets  of  functions.   Experience  with  problems  of  the  general 
type  considered  in  this  thesis  suggest  that  only  certain  particular 
choices  for  the  expansion  functions  and  for  the  testing  functions  yield 
equations  which  are  amenable  to  solution.   In  particular,  when  the  test- 
ing functions  are  chosen  to  be  Dirac  delta  functions,  the  resulting 
technique  is  usually  referred  to  as  point  matching. 

B.   EXPANSION  FUNCTIONS  FOR  THE  UNKNOWN  CURRENT 

There  are  several  popular  choices  for  the  expansion  functions  which 
have  been  considered  in  the  literature.   Some  of  the  more  common  of  these 
are  piecewise  constant,  piecewise  linear,  piecewise  parabolic  and  piece- 
wise  sinusoidal  [Refs.  20,  21].   All  of  the  aforementioned  fall  under  the 
general  label  of  the  method  of  subsections.   This  method  has  the  antenna 
contour  divided  into  a  number  of  intervals  or  subsections  where  the 
boundary  condition  of  equation  (1)  and  (2)  is  to  be  applied.   The  lengths 
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of  the  intervals  need  not  all  be  equal.   For  certain  types  of  problems 
it  could  be  advantageous  to  weight  the  lengths  of  the  intervals.   In  order 
to  keep  the  formulation  as  simple  as  possible,  however,  the  present 
solution  uses  N-l  intervals  of  equal  lengths  and  two  intervals  of  half 
the  length  of  the  N-l  intervals.   These  half  subsections  were  to  allow 
a  match  point  to  be  placed  at  the  ends  of  the  antenna  and  one  on  the 
bend  of  the  antenna.   In  addition  to  these  piecewise  functions  is  another 
family  of  expansion  functions  which  are  commonly  called  entire-domain 
expansions.   A  discussion  of  the  use  of  these  entire-domain  functions 
will  not  be  attempted  in  this  thesis. 

All  of  the  calculations  and  results  presented  in  this  thesis  are 
obtained  by  the  method  of  subsections  in  conjunction  with  point  match- 
ing.  In  particular  a  series  of  piecewise  constant  and  piecewise  sinu- 
soids are  used  to  represent  the  unknown  currents.   The  resulting 
equations  are  forced  to  be  satisfied  at  a  set  of  discrete  points  called 
match  points.   Using  Dirac  delta  functions  as  weighting  functions  ef- 
fectively requires  that  equation  (24)  be  satisfied  at  N  points  on  the 
contour  of  the  antenna. 

Harrington  and  Mautz  [Ref.  2]  have  solved  the  single  straight  wire 
problem  using  three  different  procedures,  all  within  the  scope  of  the 
method  of  moments.   In  order  of  increasing  computational  complexity, 
these  include: 

(1)  Point  matching  with  pulses  as  basis  functions 

(2)  Point  matching  with  piecewise  linear  basis  functions 

(3)  Galerkin's  procedure  with  piecewise  linear  basis  functions. 
As  the  number  of  subsections  increases,  the  second  solution  converges 
about  twice  as  fast  as  the  first.   With  antenna  segments  less  than  X/10 
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in  length,  no  significant  difference  in  results  is  observed  between 
Galerkin's  procedure  and  the  second  procedure.   Using  piecewise  sinu- 
soids is  equivalent  to  the  second  method  in  speed  of  convergance  and 
and  accuracy,  but  piecewise  sinusoids  do  not  require  the  segment  length 
to  be  less  than  A/10  in  length  to  obtain  this  accuracy  and  speed. 
Using  Dirac  delta  functions  as  weighting  functions  gives 

W  =  6(z  -  z  )U  (33a) 

m  m  z 

W  =  6(s  -  s  )U  (33b) 

s  m  s 

in  which  U  is  a  unit  vector  directed  along  the  positive  z-axis  and  U 
is  a  unit  vector  directed  along  the  positive  s-axis.   Using  the  piece- 
wise  constant  subsection  gives  the  following  expressions  for  current 

N 

I(z)  =J  IP  (z)U  (34a) 

L.    n  n    z 
n=l 

I(s)  -  T  I  P  (s)U  (34b) 

i  q  q   s 

q=l  J 

where  1   and  I   are  the  complex  constants  to  be  determined  and  the  pulse 
n      q 

functions  are  constant  over  the  increments.   This  statement  in  equation 

form  is 

1   z  <z<z   . 
Pn(z)  =  I  n    n+1  (35a) 

0  elsewhere 

1  s  <s<s   - 
P  (s)  =   <      q    q+1  (35b) 

0   elsewhere 

Similarly,  the  piecewise  sinusoid  expression  gave  the  following  current 

equations 

N  I   sin  k(z'   -z')  +  I  in  sin  k(z'-z') 
I(z)  =  I   —2 2±i 2±1 5_  Pn(z)Uz    (36a) 

n=l  sin  k(z' , ,-z1 ) 

n+1  n 
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Q   I   sin  k(s'   -s')  +  I ...  sin  k(s'-s') 

Kb)  =  I   -3 q-±± -3±i 9—  P  (s)U  .   (36b) 

=  1  Q    s 

q  sin  k(s  '  ,,-s') 

q+1   q 


C.   LOCATION  OF  MATCH  POINTS 

There  has  been  little  quantitative  study  done  on  the  precise  effect 
of  the  location  of  the  match  points  on  the  solution;  however,  it  has 
been  the  experience  of  this  author  that  the  location  of  these  match 
points  significantly  affects  the  accuracy  and  convergence  of  the  solution. 
The  number  of  match  points  and  their  location  depends  on  such  factors  as 
the  form  of  the  integral  equations  (Pockington  or  Hallen) ,  the  type  of 
current  expansion  used  (pulse,  sinusoids)  and  the  allowable  time  and 
core  for  computation. 

Extensive  calculations  have  been  carried  out  [Ref.  22]  where  the 
match  points  are  located  at  ends  of  each  current  subsection.   The  results 
of  these  calculations  indicate  the  solutions  are  quite  stable  and  reason- 
ably accurate  even  for  problems  where  a  large  number  of  match  points  is 
used. 

Other  investigations  have  indicated  that  a  good  choice  for  the  loca- 
tion of  the  match  points  is  at  the  center  of  a  current  subsection.   As 
the  author  discovered  these  conclusions  have  to  be  tempered  with  other 
information.   One  question  is  how  and  what  kind  of  experimental  verifi- 
cation is  used.   If  a  scattering  cross  section  is  used,  the  demonstrated 
results  may  differ  from  those  obtained  for  an  excited  antenna.   The  reason 
is  the  cross  section  calculations  effectively  average  the  source  and 
current  distributions  on  the  structure.   Thus,  a  convergence  study 
pertaining  to  input  admittance,  using  scattered  results,  may  be  somewhat 
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inconclusive  since  the  segmentation  and  match  point  dependence  are 
directly  related  to  source  model  and  numerical  convergence. 

All  of  the  results  presented  in  this  thesis  are  the  result  of  locating 
the  match  points  at  the  center  of  a  current  subsection  for  piecewise  con- 
stant expansion  and  at  the  ends  of  the  current  subsections  for  piecewise 
sinusoid  expansion. 
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V.   APPLICATION  OF  THIN  WIRE  THEORY  TO  AN 
ARBITRARILY  BENT  WIRE  MONOPOLE 

A.  GENERAL  CONSIDERATIONS 

A  general  theory  for  determining  the  currents  induced  on  an  arbitrar- 
ily bent  wire  monopole  has  been  outlined.   The  purpose  of  this  section  is 
to  illustrate  how  to  obtain  the  exact  solution.   The  configuration  to  be 
studied  is  composed  of  two  straight  nonorthogonal  elements  intersecting 
at  a  single  junction  which  are  driven  above  a  perfect  infinite  ground 
plane.   The  problem  is  developed  in  a  manner  which  allows  general  var- 
iations in  many  of  the  fundamental  parameters.   For  example  the  lengths 
of  the  individual  elements  may  vary,  the  angle  between  the  individual 
elements  may  vary  and  the  radii  of  the  elements  may  vary.   It  should  be 
noted  that  all  variations  in  parameters  are  still  subject  to  the  basic 
assumptions  and  restrictions  imposed  by  thin  wire  theory.   These  restric- 
tions do  not  greatly  inhibit  the  usefulness  of  the  investigation  since 
much  general  knowledge  can  be  obtained  from  the  thin  wire  approach. 

B.  COUPLED  INTEGRAL  EQUATIONS 

The  coupled  integral  equations  for  magnetic  potential  which  were 
developed  in  a  previous  section  are 


ik 
A  (z)  =  a   cos  kz  +  b   sin  kz  +  —    E  (n)  sin  k(z-n)dn 
z       z  z  to  J   z 

n=z 
z  o 

1  (  d 

[  -r1  A  (n)]  cos  k(z-n)dn  (37a) 

dy  y 

n=z 
o 
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ik 

A  (s)  =  a   cos  ks  +  b   sin  ks  +  - —  I  E  (n)  sin  k(s-n)dn 
s       s  s  w  J   s 

n=s 


-  r  I  f^U  (n)]}  cos  k(s-n)dn 

k  j  l3w  w    J 


(37b) 


n=s 


and  the  scalar  potential  equations  are 


-  J-~  d>    (z)  =  -ka   sin  kz  +  kb   cos  kz  +  k    It- [A  (n)  ]  [  sin  k(z-n)dn  (38a) 

w  z  z  i  dy  y 


n=z 


and 


'k2 

- —  6    (z)  =  -ka   sin  ks  +  kb   cos  ks  +  k 
w  s  s 


J 

n=s 


{-~[A   (n)]}  sin  k(s-n)dn.  (38b) 


3w  w 


The  constants  a  ,  b  ,  a  ,  and  b  must  assume  values  which  render  equa- 
z   z   s      s 

tions  (37)  and  (38)  in  compliance  with  the  following  boundary  conditions: 

(1)  The  scalar  potential, 4> ,  must  be  continuous  along  the  wire 
segments  except  at  the  delta  gap  generator  and  in  particular  the  scalar 
potential  must  be  continuous  at  the  wire  junction. 

(2)  The  magnetic  potentials,  A  and  A  ,  must  be  continuous  along 
the  wire  segments  and  at  the  wire  junction. 

(3)  Current  must  vanish  at  the  free  end  of  the  wire. 

(4)  The  axial  currents  must  satisfy  Kirchoff's  current  law  at 
the  wire  junction. 

The  magnetic  potentials,  A  and  A  ,  developed  earlier,  are, 


z      4tt 


-jkR 


(39a) 


and 


M3)   "  A7 
S         HTT 


-jkP 


(39b) 
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where  R  =  |r  -  r' |  and  P  =  |p  -  p' | .    R  and  P  were  evaluated  at  r'  and 
p'  equal  to  the  radius  of  the  respective  wire  segments  and  <j>'  equal  to 
zero.   Using  these  values  for  |r  -  r' |  and  J  p  —  p  *  [  and  assuming  the  sur- 
face current  flows  as  an  equivalent  line  current  along  the  s  and  z  axes 
is  equivalent  to  using  the  reduced  kernel.   This  will  be  proven  in  Ap- 
pendix A  along  with  conditions  when  the  approximation  is  no  longer  valid. 

Looking  back  at  the  development  of  equation  (37)  it  is  seen  that  the 
mathematical  formulation  was  given  without  any  physical  insight  as  to 
where  in  space  the  equation  was  valid.   In  going  from  equations  (10) 
and  (11)  to  equation  (18)  there  is  a  built-in  assumption  not  expressly 
stated.   This  assumption  is  that  the  radius  vector  is  fixed  at  some  value 
which  results  in  equation  (37)  being  valid  only  along  a  line  parallel  to 
z  and  s  directed  wires.   This  result  will  be  fully  shown  in  Appendix  B. 
Since  equations  (37)  are  valid  along  any  line  parallel  to  the  z  and  s 
directed  segments,  one  can  require  them  to  be  along  the  surface  of  the 
wire  segments  which  is  along  a  line  (a,o,z)  and  (a,o,s)  respectively. 
This  surface  is  also  where  the  tangential  E  field  is  zero.   Since  the 
source  points  and  field  points  could  not  lie  along  the  same  line,  the 
reduced  kernel  has  to  be  used  with  its  equivalent  line  current  flowing 
along  the  axis  of  the  two  wire  segments. 

Figure  1  defined  the  co-ordinate  system  which  is  composed  of  two 

co-ordinate  subsystems.   They  are  the  (X,Y,Z)  and  (P,W,S)  co-ordinate 

systems.   The  kernels  for  the  wire  segments  are  then  expressed  in  these 

co-ordinate  systems. 

-jkRZ  .  -jkRZ 

K    (x,y,z,z')    =  £— —  K^(x,y,z,z:)    =  ^— — -  (40a) 

z  R  RT 

i 

-jkRS  .  -jkR? 

K    (x,y,z,s')    =  £— —  K^(x,y,z,s!)    -  ^— —  (40b) 

S  RS  S  X  R 

l 
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A 

Us  = 

A 
A 

(4  = 


A  A 


Z 

I 

<f  jit 


0. 


Cos  9     -hUy  Sia/  & 
A  A 

us  S/v  e    +•  Uw  Cos  © 
A  A 

lie  Co j  £    —  Uw  S/'/v  d 


USjl  -  -  U3  Cos**©  +-  Uw  Si  a/  2  * 


Figure  2.   Antenna  Coordinate  Systems 
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•1      ' 
-ikr 

G    (p,w,s,z')    -  -: — — 


Gs(p,w,s,s') 


-jkr^ 


•i     z 
-jkr. 

^i /                  in        ei 
Gz(p,w,s,Zp    = — 

r. 

l 

•i     s 
-ikr. 

G    (p,w,s,s!)    = 

s    r  *    '    '    1  s 


(40c) 


(40d) 


.n2 


r-  {(x  -  x'T  +  (y  -  y')    +  (z  -  ?:r\ 


»n2 


1 


2  2  2    2 

r  =    { (p   -  p')      +    (w  -  w1)      +    (s   -   s')    }      ,    in   general. 

In   order  to   account    for  the   images    due   to   a  perfect   ground  plane   one 
must   include   the   image    current  kernels  which   are   denoted   above  by   the 
superscript   i. 

With   reference   to   figure   2, 

hz  ° 

A(7)    =   U    {[    I    (z')K   (x,y,z,zf)dz'   +    I    I    (z!  )K1(x,y ,  z,  z!  )dz!  } 

zlJz  Z  JZ1Z  11J 

z'=0  z!=-h 

h  nl  Z 

f    S  f 

+  U  I    (s')K   (x,y,z,s')ds'    -   U    .         I    (s  !  )KX(xsy ,  z,s  !  )ds  '.       (41) 

Sjo  S  SlloJLS  11 

s'=0  s!   =-h 

i  s 

and 


A(p)    =  U    {     I    (z')G    (p,w,s,z')dz'    +        I    (Z!)G1(p,w,s,z!)dz:} 

£a  I  £j  Zj  £j  JL  C*  J*  ,1 


h 

r    s 


z'=0 


z]=-h. 


+  U      II    (s')G    (p.w.s.s^ds'    -   U    .         I    (s'.)G1(p,w,s,s!)ds!.    (42) 


s  J      s 
s'=0 


SI 


SIS 


s!=-h 

1         s 


When  the  U  .  and  U  unit  vectors  are  broken  into  their  respective  U 
si      s  z 

and  U  components,  they  give  the  required  A  (r)  and  A  (r) .   Similarly 
y  z        y 

breaking  U  and  U  .  into  their  respective  U  and  U  components  gives 

Z        SI  s        w 

A  (p)  and  A  (p).   These  four  magnetic  potentiaDs  are  needed  in  equation 
37.   These  four  required  magnetic  potentials  are: 
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41T    A    /   \ 

—  A  (r) 
V      z 


h 

r  s 


I  (z')k  (x^z.z^dz1  +  cos9 


Is(s')ks(x,y,z,s')ds\  (43a) 


z'=0 


s'=0 


4^    A    (        \ 

—  A  (r) 

y  y 


sine    I  (s')k-7  (x,y,z,s')ds 
s'=0 


(43b) 


—  A  (p) 
y   s 


I  (s'){G  (p,w,s,s')ds'  +  cos  6  GX(p,w,s,sI)}ds'  + 
h 

;   z 


cos  e 


I  (z')G  (p,w,s,z')dz' 


and 
4tt 


z'=0 
h 


A  (p) 

y   w 


=  -sin  26 


Is(s')Gs(p,w,s,s')ds 


s'=0 


(43c) 


sin  6    I  (z')G  (p,w,s,z')dz' 
z'=0 


(43d) 


Placing  equations  (43)  into  equations  (37)  gives 


h  h 

(  z  t   s 

t^-     I  (z')KZ(x,y,z,z,)dz'  +  cos  6    I  (s  '  )RS  (x,y ,  z,  s '  )ds '    -  a  cos  kz  + 

't  It  ^  fa  t* 


z'=0 


b  sin  kz  +  — 

Z  CO 

and 
h 


s'=0 
z 


E  (n)sin  k(z-n)dn  -  — 
z  K 


{-^-[A  (n)]}  cos  k(z-n)dn    (44a) 

dy   y 


n=z 


n=z 


i 
I  (s'){G  (p,w,s,s')  +  cos  G  (p ,w,s ,s ' )}ds'  +  cos( 


h 
-  z 


I  (z')G  (p,w,s,z')dz 


s'=0 


z'=0 


jk 


a  cos  ks  +  b  sin  ks  + 
s  s  0) 


E  (n)sin  k(s-n)dn  -  -   {-A  A  (n)]}cos  k(s-n)dn. 

s  K  j  l3w  w 


n=s 


n=s 


(44b) 


Examination  of  equations  (44)  indicated  that  the  partial  derivatives 

of  A  (r)  and  A  ( p)  are  needed.   The  results  for  A  (r)  will  be  shown  and 
y        w  K  y 

A  (p)  will  follow  in  a  similar  manner, 
w 
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The  form  of  A  (r)  is 


and 


-jkR 
Ay(r)  -  £-_  (45a) 

3A  (r)    3A  (r)  3R 
y    =   y   

3y       3R   3y 

9Ay(r)     (1  +  jkR)(y  -  y')e"jkR  .  (45b) 

9y   "  R3 


Similarly,  for  A  (p  )  , 


w 


3A  (p )     ,..    ,    .,  v  ,  ,N  -jkr  ..,. 

w    _    (1  +  jkr) (w  -  w' )e  J  (46) 

3w  3 

r 

Placing  equations  (45)  and  (46)  into  equations  (44)  and  examining 

the  results  closely,  it  is  seen  that  the  double  integrals  have  an  inner 

integral  of  the  form 

[^|  +  ~]e~jkR  cos  k(B-n)dn  (47) 

J   R    R 

1 

where       R  =  {a2  +  g(a)  +  [n  -  f(a)]2}2 

and  f  and  g  are  any  function  of  the  parameter  a.  Reference  [23]  gave 
the  exact  solution  for  this  integral  which  reduced  the  double  integral 
to  a  single  integral  and  thereby  greatly  reducing  the  computer  inte- 
gration time  and  numerical  errors.  The  exact  derivation  will  be  shown 
in  Appendix  C.  Now  place  the  result  of  Ref.  [23]  into  equation  (44a). 
With  rearrangement  of  equation  (44a)  the  following  equation  is  obtained; 
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[— ]    a     cos   kz  +    [— ]   b      sin  kz  +    [-^J^] 
y        z  y        z  y  w  z 


n=z 


4tt 

yk 


h 

;  s 


, .  r [n-f (s ' ) ] cos  k(z-n)    -   jRsin  k(z-n) 


I    (s')     (Ky(x,y>nts'){Ln-US,;jCO%K^-n;    -   JKS"n   KKZ~n)}      )ds» 

,    J  R(RZ   -    [n-f(s')]Z) 

s   — s  n=z 

h  h 

z  r   s 

I    (z')KZ(x,y,z,z')dz'    -   cos    6        I    (s ' )KS (x,y , z,s ' )ds '    =  0  (48) 

z  J       s 

z'=z  s'=s 


Similarly,    for  equation    (44b): 


[ — ]    a      cos   ks  +    [ — ]   b      sin  ks  +    [ — — ] 

y        s  '    y        s  yw 


E    (n)    sin  k(s   -  n)dn  + 
s 


n=s 


4-rrsin   20    (        .    M    />i/     T7  „   „  ,  A  j  [n-f  (s  '  )  ]  cos   k(s-n)    -   jrsin  k(s-n)i 


s 


li 


S'=S 


r(r2   -    [n-f(s')]2) 


)  ds  '    + 
n=s 


4irsin  9     {  i  (S»)   (G^P,w,n,z'){^-^2,)]co%k(s--n)  "  Jrs^-Miiil>}    )d 


r(r2  -  [  -f(z')]2) 


s  '  - 


n=s 


h 

r   s 


I  (s')[G  (p,w,s,s')  -  cos  G  (p,w,s,s ' ) ] ds 
s       s  s 


s'=s 


cos 


I  (z»)G  (p,w,s,z')dz'  -  0   . 


(49) 


For  a  structure  driven  by  a  delta  gap  source  of  emf,  V,  on  an  element 
at  a  location  displaced  g  units  above  a  ground  plane  the  following  equa- 
tions are  obtained  [Ref.  20]: 


V  (z)  =  &■   I  -Vfi(n-g)  sin  k(z-n)dn 


(50a) 


n=z 


Vz(z) 


_  Ji 


[V  sin  k(z-g)  U(z-g)]. 


(50b) 
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It  is  desirable  to  obtain  a  representation  of  the  emf  when  the 
structure  is  base  driven,  i.e.,  when  g  is  zero.   A  structure  above  a 
ground  plane,  regardless  of  the  nature  of  its  excitation,  always  gave 
rise  to  zero  scalar  potential  everywhere  on  the  ground  plane.   This  is 
a  direct  result  of  the  uniqueness  theorem  and  boundary  conditions  at  a 
perfect  conducting  surface.   Thus,  it  is  necessary  that,  if  the  element 
and  its  source  are  to  be  equivalent,  the  impressed  emf  must  be  sym- 
metric; otherwise,  the  currents  and  charges  in  a  free-space  "equivalent" 
of  the  structure  above  a  ground  would  not  satisfy  imagy  theory  require- 
ments . 

It  is  from  this  discussion  that  the  modification  to  equation  (50b) 
is  found.   When  a  base-driven  structure  is  being  analyzed,  one  sets  g 
equal  to  zero  and  V  equal  to  V/2.   If  V/2  were  not  used  in  place  of  V, 
one  would  be  treating  an  antenna  driven  by  2V  when  the  free-space 
"equivalent"  structure  is  formed.   These  results  are 

V  (z)  =  -  Ji  V  s.n  kz  ^  (51) 
z         to  2 

Since  the  s-directed  segment  is  not  driven  by  a  separate  generator,  this 

gives 

s 

V  (s)  -  0  or  Jii   [0]  sin  k(s  -  n)dn  =  0  .  (52) 


S  CO 


n=s 
o 


Now  the  boundary  conditions  are  applied.  Due  to  symmetry  of  the 

charges  on  the  structure  and  its  image,  the  scalar  potential  must  be 

zero  at  the  point  where  the  vertical  element  enters  the  ground  plane. 

This  implies 

.,2 

-  ^4>Z(0)  =  0  (53a) 

which  gives  b   =  0  (53b) 
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Next   by   condition  one      the  scalar  potential  must   be   continuous   at   the 
wire  junction  which   gives 


<j>z(h  J  =  4>s(0) 


(54a) 


and  results  in 


-ka  sin  kh  -  kb  +  k   {— —  [A  (n)]}sin  k(z-n)dn  - 

z    z    s    J  3y  y 

n=z 


J^  ^  cos  kh  =0  (54b) 
a)  2       z 


Conditions  three  and  four  require 


1=0   and 


(55) 


IZ  -  1S  . 
N     1 


(56) 


Clearing   all   constants   and  applying  all  boundary   conditions,    the   fol- 
lowing equations    for   I    (z)    and   I    (s)    are   obtained: 

z  s 


h 

U  f  s 

[— ]a   cos   kz  -        I    (s')(Ky(x,y,n,s'){ 

s'=0 
h 

•17  f     z 


,  ^ ;  [n-f(s')]cos   k(z-n)    -  jRsin  k(z-n) 
R(R2   -    [n-f(s')]2) 


)ds' 
n=0 


h 

r   s 


sin  kz  -        I    (z')K   (x,y ,z, z' )dz'    -   cos 
z'=0 


I    (s')K   (x,y,z,s')ds 


s»=0 


-   0 


(57a) 


[ — 1  a   cos   ks   + 

y     z 


h 
r  s 


, > r [n-f (s ' ) ]cos  k(s-n)    -   jrsin  k(s-n) 


r(r2   -    [n-f(s')]2) 


}hdS' 


.471-  r  z 


I8(B,)(G1(p,w,n,B,){- 

h,  s 

z,  ,  N  r  [n-f  (z'  )  ]cos  k(s-n)    -   jrsin  k(s-n) 


-K-^-]b   sin  ks  +        I    (z')(GZ(p,w,n,z')l 

[,:Q  r(r2-    [n-f(z')]2) 


}   )dz' 


h 
f  s 


I    (s')G    (p,w,s,s')ds'    -   cos 


s'=0 


h 

'   s  i 

I    (s')G    (p,w,s,s')ds 
s  s 

s'=0 


-   cos 


I    (z')G   (p,w,s,z')dz'      =     0 


(57b) 


z'=0 
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and 


h  h 

s 


f     S  7 

r47rn         .      ,,        ,        _     .    ,*  rvy .  ,  N  r  [n-f  (s  '  )  ]cos   k(h   -n)-jRsin  k(h  -n)  i  ^  , 

-    [ — ]a  sin  kh     +        1    (s')lK   (x,y,n,s'){ z — '—^ i-« — '-\  d 

u      z  z  s  „  „  '  J 


s'=0  R(R2   -    [n-f(s')]2)  n=0 


-    [— ]b      -  -^  cos   kh        =     0  (57c) 

y      s        60  z 


Iq     -     0  (57d) 

XN      =      A    '  (57G) 

There   are    five   unknowns    a    ,    b    ,    a    ,    I    (z1),    and   I    (s')    to  be    found 

s        s        z        z  s 

and  five  equations  which  relate  these  quantities;  hence,  there  is  a  unique 
solution  to  this  system  of  equations  for  the  five  unknowns .  [Ref.  20]. 
Obtaining  the  current  distributions,  all  the  other  desired  quantities 
such  as  electric  field,  input  impedance  can  be  determined. 

C.   NUMERICAL  SOLUTION 

The  system  of  integral  equations  given  in  the  previous  section  must 
be  solved  in  order  to  determine  the  unknown  current  distributions.   The 
direct  integration  technique  as  outlined  earlier  is  used  to  effect  a  nu- 
merical (discrete)  solution  for  the  current  distributions. 

The  first  step  in  the  direct  integration  technique  is  the  choice  of 
the  expansion  functions  and  the  testing  functions.   As  pointed  out  in  an 
earlier  section,  the  method  of  subsections  in  conjunction  with  point 
matching  is  a  very  convenient  approach.   The  testing  functions  are  Dirac 
delta  functions  since  point  matching  is  used. 

The  expansion  functions  for  the  unknown  currents  are  chosen  to  be 
piecewise  constant  for  its  ease  of  use  and  piecewise  sinusoids  for  their 
improvement  in  convergence  and  required  computation  time.   Symbolically, 
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the  piecewise  constant  function  is  represented  by 


N 

I  (z')  =  I    IZ[U(z'  -  z')  -  U(z'  -  z'   )]  (58a) 

z        -,  n         n  n+1 

n=l 

I  (s»)  =  V  IS[U(s'  -  s')  -  U(s'  -  s'   )]  (58b) 

S        *   q  q  q+1 

and  where  U(a)  is  the  familiar  unit  step  function 


1   a  >  0 
U(a)  =  <  (59) 

0   a  <  0   . 


Consistent  with  the  above  description,  the  current  can  be  written  as 

N-l 

I  (z')  =  IZU(z!  -  z')  +  I   IZ[U(z'-z')  -  U(z'-z'   )]  +  IZU(z'-z')    (60a) 
z        1   I  n  n       n  n+1      N      N 

n=2 

and 

Q-l 

I  (s«)  =  IZU(s!  -  s«)  +  I   IS[U(s'-s')  -  U(s'-s'   )]  +  I^U(s'-s')  .  (60b) 
s         12  ^2  q        q  q+1       Q       Q 

Similarly,  one  can  write  for  the  piecewise  sinusoids  the  following  ex- 
pressions: 

N  IZsin  k(z'-z')  +  IZ   sin  k(z'-z') 

i  (z»)  =  I  -2 °±i 2+i ^-{u(z'-z;)  -  u(z'-z;+1)} 

n=l  sin  k(z' , ,-z1) 

n+1  n  rr-\    \ 

(61a) 

Q  ISsin  k(s'   -s')  +  IS,.,sin  k(s'-s') 
I  (s')  -  I   -3 3+1 3+1 !_{u(s._s.)  _  u(s._s.   )} 

S      q=l  sin  k(s'   -s')  q  q+± 

q+1   q 

(61b) 
With  the  assumption  that  the  segments  were  all  equally  spaced,  the  term 

z'-zf  is  a  constant.   This  allows  a  term  delta  (A)  to  be  defined  as 
n+1  n 

A  =   ,    1 .  (62) 

sin  k  (constant; 
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With  this  definition  and  equations  (61)  the  current  can  be  written  as 

1Z  IZ 

I  (2)  =  -i  sin  k(z'-Z')U(Z'-zT)  +  -^tl  sin  k(z  ' -z' )U(z '-z  ' )  +      (63) 

z        A         LI  A  n        N 

1  N 

f-  7  IZ{sin  k(z'-z'  jflKz'-z'  1)-U(z,-z*))+  sin  k(z'  ,  -z')  (u(z'-z' )-U(z f-z'  j)} 
A  ^_  nl  n-1  *"      n-1        n  '  n+1     v      n        n+1  J  J 

n=2 


and  similarly  for  the  s  current.   With  these  expansion  functions,  the  cur- 
rent distribution  appeared  graphically  as  shown  in  Figures  3-9. 
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BajaJiuslTTlW 


Figure  3.   Real  and  imaginary  parts  of  the  current  distributions  for  the 
radius  =  .001588X,  theta  =0°,  36  segments. 
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Figure  4.   Real  and  imaginary  parts  of  the  current  distributions  for  the 
radius  =  .004763A,  theta  =0°,  20  segments. 
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Figure  5.   Real  and  imaginary  parts  of  the  current  distributions  for  the 
radius  =  .001588X,  theta  =  0°,  30°,  45°,  60°,  85°,  90° ,  36  segments. 
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Figure  6.   Real  and  imaginary  parts  of  the  current  distribution  for  the 
radius  =  .003175X,  theta  =  0°,  30°,  45° ,  60°,  85°,  90°,  32  segments. 


45 


ssjadoBTtTTW 


o  K 


oajadostix™ 


Figure  7.   Real  and  imaginary  parts  of  the  current  distribution  for  the 
radius  =  .004763X,  theta  =  0°,  20°,  45°,  60°,  85°, 90°,  25  segments. 
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Figure  8.   Real  and  imaginary  parts  of  the  current  distributions  for  the 
radius  =  .001588A,  .003175A,  .004763A,  theta  =  45°. 
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Figure  9.   Real  and  imaginary  parts  of  the  current  distributions  for  the 
radius  =  .001588A,  .003175X,  .004763A,  theta  =  90°. 
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VI.   RESULTS 

The  system  of  integral  equations  and  the  proposed  solution  techniques 
are  quite  general,  and  they  allow  variations  in  many  of  the  parameters  in 
the  basic  configuration.   However,  certain  parameter  variations  appear  to 
be  more  interesting  than  others.   This  study  concentrates  its  efforts  in 
three  main  areas: 

1.  The  ability  of  the  numerical  solution  to  give  the  result  obtained 
by  King  in  Ref.  26  when  the  angle  between  the  two  wires  is  zero. 

2.  Current  as  a  function  of  the  angle  between  the  two  wires. 

3.  Current  as  a  function  of  wire  diameter. 

The  antenna  is  assumed  to  be  driven  by  a  delta  gap  generator.   The 
current  distributions  are  presented  in  terms  of  a  dipole  vice  a  monopole. 
The  real  and  imaginary  parts  of  the  current  are  plotted  for  various  bends 
of  the  antenna.   The  total  length  is  A/4  for  all  the  figures  with  the  an- 
tenna bent  at  A/8  to  give  two  equal  lengths.   Although  these  results  are 
the  only  ones  presented,  the  computer  program  allows  the  antenna  to  be 
bent  at  any  point  and  the  antenna  to  have  any  length. 

Figures  3  and  4  show  the  current  for  two  different  radii  at  an  angle 
of  zero  degrees  for  a  comparison  against  King's  frill  feed  in  Ref.  26. 
These  figures  also  show  the  result  from  Ref.  27  and  28  against  King. 

References  27  by  D.C.  Kuo  and  B.J.  Strait  at  Syracuse  University  uses 
a  modified  Pockington  formulation  with  thin-wire  approximations.   4>  is 
used  instead  of  A  and  derivatives  are  approximated  by  finite  differences 
over  the  same  intervals  used  for  integration.   Galerkin's  procedure  is 
used  with  triangle  expansion  functions  which  are  non-zero  only  over  four 
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consecutive  segments.   Successive  triangles  overlap  every  two  segments 
except  on  the  ends  of  open  wires.   Junctions  are  treated  as  open  wires 
each  having  a  two  segment  overlap  at  the  junction.   An  extended  operator 
is  used  because  the  first  order  derivatives  of  the  triangle  functions 
do  not  exist  at  certain  points. 

Reference  28  by  MB  Associates  and  called  the  Antenna  Modelling  Pro- 
gram uses  a  Pockington  formulation  with  thin-wire  approximations.   Point 
matching  is  used  with  the  current  being  expanded  as  a  sine,  cosine  and 
constant  with  an  interpolation  procedure  so  only  N  linear  equations  need 
to  be  solved.   The  extrapolated  current  for  a  given  segment  is  forced  to 
match  the  center  current  values  in  two  adjacent  segments.   This  interpol- 
ation procedure  is  also  applied  to  wire  junctions. 

Figures  5,6   and  7  give  the  currents  for  antennas  of  three  different 
radii  having  0,  30,  45,  60  and  90  degree  bends.   King's  input  values  for 
a  0  degree  bend  are  also  plotted.   Figures  8  and  9  show  the  currents  for 
different  radii  antennas  plotted  on  the  same  graph  for  bends  of  45  and  90 
degrees,  respectively. 

Several  general  comments  concerning  the  data  presented  may  prove  help- 
ful to  the  reader.   All  the  results  presented  here  are  for  the  case  of  a 
one  volt  drive  and  assumes  an  open-ended  perfectly  conducting  hollow 
cylindrical  tube  as  the  antenna.   Extensive  calculations  have  been  carried 
out  for  bends  of  0  to  90  degrees  for  each  antenna  radius  given,  but  the 
plotting  of  all  these  angles  on  a  single  graph  would  obscure  the  overall 
view.   It  is  believed  that  the  data  given  is  sufficient  to  give  the  reader 
an  indication  of  the  general  variation  of  the  current  with  respect  to  the 
parameter  of  interest. 
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VII.   CONCLUSIONS 

In  this  thesis  a  new  method  is  presented  for  analysis  of  an  arbitrarily 
bent  monopole.   The  arbitrarily  bent  monopole  is  a  good  choice  because  it 
contains  cross-  coupling  terms  in  the  Hallen  formulation  which  are  not 
present  in  the  0  or  90  degree  case.   The  currents  had  to  be  obtained 
using  the  appropriate  boundary  value  problem  based  on  Maxwell's  equations 
which  are  formulated  as  a  operator  equation.   The  general  procedure, 
"method  of  moments",  is  used  to  represent  the  operator  by  a  matrix.   Geo- 
metrical symmetries  are  utilized  to  simplify  the  matrix  inversion. 

Numerical  results  are  compared  with  data  obtained  by  different  theo- 
retical techniques  and  with  experimental  results  where  possible.   Agree- 
ment is  generally  very  good  and  this  is  interpreted  as  an  indication  of 
the  accuracy  of  the  solution. 

Certain  parameters  should  now  be  concentrated  on  if  the  conclusions 
are  to  be  used  properly.   These  efforts  should  be  in  these  areas: 

1.  Choice  of  bases 

2.  Antenna  drive 

3.  Discontinuity  at  the  bend. 

As  can  be  seen  by  equation  (36)  the  sinusoid  current  expansion 
equations  are  continuous  but  their  derivatives  are  not.   This  implies 
mathematically  that  there  are  pockets  of  charge  at  the  end  of  each 
subinterval  which  physically  are  not  there.   If  one  is  interested  in 
the  near  field,  these  mathematical  charges  cause  discontinuities  in  the 
near  field.   One  solution  has  been  to  use  a  least  square  fit  to  the  cur- 
rent values  obtained. 
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The  delta  gap  drive  was  chosen  because  it  was  an  ideal  model.   Even 
though  the  delta  gap  gives  a  physically  meaningless  infinite  input  sus- 
ceptance  associated  with  the  knife  edges,  the  answer  is  close  to  that 
achieved  by  experiment.   The  reason  is  that  any  method  that  approximates 
the  current  by  a  series  of  continuous  functions  will  not  reproduce  the 
extremely  sharp  rise  in  amplitude  near  z  =  z'  unless  a  very  large  number 
of  terms  is  used.   A  better  drive  would  be  the  frill  feed  as  given  in 
reference  [29]  which  corresponds  closely  to  the  actual  drive  used  on 
physical  antennas. 

Although  not  plotted  on  the  figures  given,  a  discontinuity  in  cur- 
rent was  obtained  within  one  subinterval  either  side  of  the  bend.   This 
was  to  be  somewhat  expected  because  the  formulation  considers  two  straight 
wires  joined  at  a  point  which  gives  a  sharp  bend  where  the  cross  section 
is  abruptly  changed.   When  one  uses  the  approximate  kernel,  I(z)  must  be 
a  slowly  varying  function  in  the  vicinity  of  z  =  z'  for  the  approximate 
kernel  to  be  within  an  order  of  the  square  of  the  radius  of  the  exact 
kernel.   This  is  not  true  at  the  bend.   It  is  this  slow  variation 
or  fast  variation  that  causes  the  various  methods  of  iteration  (Hallen, 
Gray , King-Middleton)  to  require  more  or  fewer  terms  in  representing  the 
current. 
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APPENDIX  A:  REDUCED  KERNEL 


Equation  (39)  in  the  thesis  must  be  investigated  and  R  and  P  deter- 
mined in  a  more  precise  manner.   A  detailed  development  of  the  vector 
potential  must  now  be  presented. 

The  vector  potential  due  to  a  volume  current  density  can  be  written 


as 


A(r) 


4tt 


ff   J(i:')e  jk'r  r'  'dV(r') 


(1) 


V 


r  -  r 


Looking  at  Figure  10,  it  is  seen  that  the  vector  r  locates  the  point  at 
which  A(r)  is  to  be  determined  and  r'  locates  the  source  point.   With  a 
symmetrically  driven  antenna  whose  wall  thickness  approaches  zero,  the 
vector  potential  has  only  a  z  component  which  can  be  written  as 


A  (r)  =  ^ 

z      4tt 


2it   a  ..  | , 

f      f  f        T  /  i   i\  -jk  r-r' 

'   '        Jz(r',z')e  J  ' 


z'  ({,'=0  r'=0 


r-r 


r'dr'd^'dz'     (2) 


or  integrating  over   r      gives 


-2tt 


A  (7)  =  -£ 

Z  4lT 


f     f 


[alz(a, z') ]e 


-jk | r-r' 


z'    (J>'=0 


r-r 


dcfi'dz 


(3) 


Now   due    to   symmetry 


or 


I(z')    =   2TraJ    (a,z') 


aJ    (a,z')    =  —  I(z')       . 

Z  /IT 


Substituting    (4b)    into  equation    (3)    gives 

2n 


A   (r)   =  -^ 

Z  4lT 


1        e 


l^tel 


-jk   r-r' 


-d* ' }  dz ' 


<J>'=0 


r-r 


(4a) 
(4b) 


(5) 
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Figure  10.   Antenna  Geometry  assuming  a  surface  current  in 
the  z  direction. 
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where  r1  is  constrained  to  remain  on  the  surface  of  the  cylinder. 

Appealing  to  basic  vector  algebra  and  trigonometry  and  being  aided 
by  Figure  10  one  can  write  r  and  r'  in  the  following  form: 

7  =  xU  +  yU  +  zU  =  r  cos  c^U  +  r  sin  <J>U  +  zU  (6a) 

x     y     z  x  y     z 

«-\  y\.  *<.  --.  /v  ^ 

r'=  x'U     +  y'U     +  z'U     =   r'cos    cj)TU     +  r'sin   <j>'U     +  z'U      .       (6b) 
x  y  z  x  y  z 

Now    | r  —   r ' |    is  written   as 

1 

|7  -   r'|    =    [(x  -    x')2   +    (y   -   y')2   +    (z   -    z')2]2 

1 

=  [r2  +  (r')2  -  2rr'cos(<j>  -  <j>»)  +  (z  -  z*)2]2      (7) 

If  |r  -  r1 I  is  evaluated  at  r'  equal  to  the  radius  of  the  z-directed 

wire  and  <j>  equal  to  zero  then 

1 
|r  -  7f  |  ,_  =  [r2  +  a2  -  2ra  cos  <j> '  +   (z  -  z')2]2  (8) 

cj>  -0 

Trigonometry  gives 

1  -  cos  2B  =  2  sin  3  (9) 

Applying  equation  (9)  to  (8)  with  r  =  a  gives 

1 
|7  -  7f  |  =  [4a2  sin2(-^)  +  (z  -  z')2]2   .  (10) 

Going  back  to  equation  (5)  and  inserting  equation  (10)  gives 
„  ,      ,   f2^-jk[4a2Sxn2(f)  +  (z  -  z')2]1/2 
z.         Ji=0   [4a  sin  (y-)  +  (z  -  z')  J 

One  sees  immediately  that  at  z  equal  z'  and  -  <f>'  equal  0  or  2tt  the  inner 
integrand  is  singular;  however,  the  integrand  is  logarithmically  singular 
and  is  integrable.   Even  though  it  is  integrable,  evaluation  of  the  in- 
tegrand by  numerical  techniques  is  extremely  difficult.   This  is  one 
major  reason  why  a  reduced,  average  or  approximate  |r  -  r'|  is  found. 
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For  |z  -  z'|  >>  a   and  for  thin  wires  where  a  <<  X    the  following 
approximation  can  be  obtained. 

|r  -  7«|  =  [(a  -  a')2  +  (2a  sin  [^])2]1/2  =  |z  -  z'  |       (12) 

Even  this  for  thin  wires  is  not  a  good  approximation  where  z  =  z'.  Ob- 
serving that  2a  sin  (— )  is  the  distance  from  a  point  (a,0,z')  to  a 
point  (a^'jz),  it  is  seen  that  both  points  are  on  the  same  circle  of 
radius  a  in  a  constant  z'  plane.   Looking  at  Figure  10  it  is  seen  that 
2a  sin  (y~)  was  a  distance  that  assumed  values  from  0  to  2a   as  <j>' 
varied  from  0  to  2tt  with   a  mean  value  of  a.   Placing  the  mean  value  into 
equation  (11)  and  integrating  gives 


2  2  1/2 

-jk[aZ  +  (z  -  z')  ]  ' 


A  (a,z)  =  -£  I(z')—  -^-        -T-T79  dz '  ■         (13> 

47r  I  [a2  +  (z  -  z')2]1/2 

It  is  restated  that  equation  (13)  is  valid  only  for  a  <<  A  and 
a  <<  h.   Now  going  back  and  using  Figure  11,  equation  (5),  and  assuming 
the  current  to  be  a  filament  along  the  z  axis,  the  vector  potential 
A  (r)  at  a  point  (a,0,z)  would  be 


z 


A  (..O,.)  -  £■ 


2  2  1/2 

-jk[az  +  (z  -  s'rr' 
I(z,) — r^— ^7172  dz?  (14) 

[a  +  (z  -  z')  ] 


which  is  exactly  the  same  equation  as  equation  (13).   This  is  a  most 
important  point.   Using  the  reduced  kernel  is  the  same  as  treating  sur- 
face current  on  the  cylinder  as  if  it  were  an  equivalent  line  current 
along  the  cylinder  axis. 
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i-r*  *au  ♦ (z - z '  )u 


(a. 0,2') 


Figure  11,   Antenna  Geometry  assuming  a  filamentary  current 
along  the  z  axis. 
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APPENDIX  B:   SOLUTION  OF  COUPLED  INTEGRAL  EQUATIONS 

Equation  (10)  and  (11)  in  the  thesis  give  the  two  coupled  partial 
differential  equation  that  must  be  solved.   Equation  (12)  gave  the  form 
that  equation  (10)  and  (11)  must  be  in  to  apply  the  solution  given  in 
equation  (17).   Equation  (10)  gives  the  form  for  the  z-directed  wire  as 

-Va2W]   +  *■%&   -  ^(r)   -   gfbjA/r)]    .  (10) 

oZ 

Notice  that  A  (r)  in  cylindrical  coordinates  depends  on  (R,<|),z)  where 
R  is  the  radius  vector  in  these  coordinates.   If  one  fixes  R  and  <J),  then 
A  (r)  will  be  a  function  of  z  alone  along  any  line  parallel  to  the  z-axis, 
Equation  (10)  then  becomes 

2  •    2 

-^[A    (z)    +  k2A    (z)    =  ^-EP(z)    -  ^[~A    (z)]     .  (1) 

2      z  z  to     z  dzdyy 

dz 
Equation    (12)    in    the   thesis    is 

d2f(x)  2  3V(x)  n9x 

Examining  equation  (12)  it  is  seen  that  equation  (1)  may  be  substituted 
directly  into  equation  (12)  and  equation  (17)  which  is  the  solution  to 
equation  (12)  may  be  applied  to  equation  (1).   Thus  it  is  seen  that  for 
equation  (10)  to  use  the  solution  given  by  equation  (17)  it  must  be  in 
equation  (1)  form.   This  means  equation  (10)  solution  is  good  only  along 
a  line  parallel  to  z-axis. 


APPENDIX  C:   REDUCTION  OF  A  DOUBLE  INTEGRAL 

In  order  to  show  the  solution  to  equation  (47)  in  general,  the  sol- 
ution will  be  merely  given  and  then  differentiated  to  show  that  it  does 
indeed  give  equation  (47)  even  with  arbitrary  signs  and  the  addition  of 
a  constant.   The  solution  is 


t      .,     ,   _  .,R  -jkR 

(Jy  +  -±)e  3      cos  k(3-e)de  = j-{  [e-f  ]  cos  (3-e)  -  jRsin  k(3-e)} 

J   R    R  R[R-(e-f)  ] 

where  (1) 


R  =  {a2  +  g2(c)  +  [e-f(C)]2J1/2. 


Now  let  the  solution  in  general  be 


-jkR 

—  {(e±f±c)cos   k(3-e).  -   jRsin  k(g-e)}  (2) 

R(R-[e±f(n)±c]    ) 


and 

R  =    {a2  +  g2(n)   +    [e±f(n)±c]2}1/2,  (3) 

where   g  and   f   are    functions    of   the   parameter  n. 
Now  taking   the   derivative   of   the   function 


-jkR 

-JL{     - — [(e±f±c)cos  k(3-e)    -  jR  sin  k(g-e)]}   = 

aZ   R  -R(e±f±c) 

d                 ~jkR 
[(e±f±c)cos  k(3-e)    -  jR  sin  k(3-e)]  -=-   [-r— - j]     + 

R  -R(e±f±c) 

-jkR 

Ht— =-]  ~   [(e±f±c)cos  k(3-e)    -  jR  sin  k(3-e)].  (4) 

R3-R(E±f±c)2      dG 
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Working   on    the   first   derivative    in   equation    (4)    gives 

~~jkR  3  2      -ikrr  +  f+c)e~jkR 

]    =    [{R  J-R(e±f±crH-^-        Cje 1   - 


dC    R3-R(e±f±c)2  R 

e-jkR  {3R2(£±f±c)    _   2R(£±f±c)    _    (e±f±c)2(E±f±c)}] 


[R3-R(e±f±c)2]2 


This    reduces    to   give   for  the    first   part   of   equation    (4) 
■jk       1  > 


rJkR(c±f±C)  P 


[- (c±f±c)      R  R      ]    [(£±f±c)    cos   k(e_e)    _  jR  sin  k(3-e)]      .         (5) 

R  -R(e±f±c) 


Similarly    for   the   second   derivative   in  equation    (4)    gives    the   following 
result 


[(e±f±c)k   sin  k(e-e)   +   cos  k(S-e)   +  jkR   cos   K(B-e)    -   J    Sln  k(g   e)  (e±f±c)  j 

K 

-jkR 

[-3— 2    ]       '  (6) 

R  -R(e±f±c) 


Adding  equations  (5)  and  (6)  and  factoring  out  the  common  elements  re- 
sults in 

-jkR 

~Y^ y[(e±f±c)k   sin  k(B-e)   +  cos  k(3-e)   +  jkR  cos   k(3-e)    -        (7) 

R  -R(e±f±c) 

2  I  1 

j(e±f±c)sin  k(B-e)        (e±f±c)    (-jk-R)cos  k(g-e)      (e±f±c)  (-jk-R)  (-jR  sin  k[g-e]) 

R  R  R 


Now  start  cancelling  common  terms 

(e±f±c)k  sin  k(g-e)  -  (e±f±c)k  sin  k(a-E)  -  0  (8a) 

j (e±f±c)sln  k(g-e)     j (e±f±c)sin  k(B-e) 


R  R 
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=  0  (8b) 


which  gives  as  a  result  after  some  algebra 
-jkR 


e 


(COS  f'"E>    {[K2-(e±ac)2][l  +  jkR]})      .  (9) 


3         2^2 
R  -R(e±f±c)      R 


Cancelling  common  factor,  one  has 

-ikR 

e  rcos  k(B-c)n  ,  i1t,»i     -jkR     /0   .  r jk  ,   1, 

I 5^ d  +  JkR)/  =  e     cos  (3-e)Hy  +  -~] 

R         R  R    R 


(10) 


which  is  equation  (47)  in  the  thesis. 

Thus  equation  (2)  and  (3)  are  required  functions  in  general  and 
equation  (1)  is  the  necessary  function  when  there  is  no  arbitrary  con- 
stant or  arbitrary  signs  in  the  integrand. 
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